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Abstract. We present an analysis of precision radial velocity measure- 
ments for 580 stars from the Keck survey. We first discuss the detection 
^ \ threshold of the survey, and then describe a Bayesian approach to con- 

strain the distribution of extrasolar planet orbital parameters using both 
detections and upper limits. 
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', 1. Introduction 

Attempts to characterize the distribution of mass, orbital radius, and eccen- 
Q ■ tricity of extrasolar planets are hampered by the lack of knowledge of detection 

\ sensitivities. Following our previous study of the Lick survey (Gumming, Marcy, 

^ ■ & Butler 1999, hereafter CMB), the aim of this work is to carefully assess the de- 

tection threshold and selection biases for a sample of stars from the Keck survey. 



and use this information to constrain the distribution of orbital parameters. 
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■ 2. Search for Companions 



We analyse radial velocity measurements for 580 F,G,K,M stars from the Keck 
survey. Typically, each star has 3 or 4 measurements per year over 2 to 5 years 
with measurement errors between 2 and 5 m/s. 

We adopt a least-squares approach. First, we check whether the velocity 
scatter is consistent with measurement errors and intrinsic stellar "jitter". The 
jitter, which arises from a combination of convective motions, magnetic activity, 
and rotation, is estimated from stellar parameters such as rotation period and 
mass. Assuming Gaussian statistics, we calculate the probability of obtaining 
as large or larger than the observed value from measurement errors and jitter 
alone, adopting a threshold probability of 10~^. The solid line in Figure 1 shows 
the signal to noise needed for a signal to be detected by this method 50% of 
the time as a function of the number of data points n. We define the signal to 
noise as K/y/2a, where K is the velocity amplitude, and a the expected noise 
amplitude. Even with only a few measurements, signals > 2a are detected. 
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Figure 1. Signal to noise ratio Kj^a needed for detection by x^- 
test (solid line) or periodogram (dashed lines). For small n, the dashed 
lines increase much more rapidly than 1 / ^Jn. 



Next, we look for periodicities using the Lomb-Scargle periodogram (Scarglc 
1982; Walker et al. 1995; CMB), which is similar to the F-statistic, measuring 
the improvement in between a model which includes a sinusoid and one that 
docs not. We look for the maximum periodogram power between 2 days and 
10 years, and use Monte Carlo methods to assess its significance. The dashed 
lines (analytic result) and squares (numerical result) in Figure 1 show the signal 
to noise ratio needed to identify a periodic signal 50% (filled squares) and 90% 
(open squares) of the time (significance level 10~^) as a function of number of 
data points. To obtain analytic estimates, we follow the approach of Horne & 
Baliunas (1986), but using the correct statistical distribution for periodogram 
powers in the presence of noise (CMB Appendix B). We find the signal to noise 
ratio needed for detection 50% of the time is Kj^fla = [(M/F)2/("-3) - 1]V2^ 
where M the number of independent frequencies searched, and F the false alarm 
probability needed for detection. For large n the 50% threshold is K/y/2a ~ 
[21n(M/F)/n]^/2, or « b/^Jn for typical values F = lO'^ and M = 1000. 

Figure 1 is striking because when looking for periodic signals, we are used 
to being able to "dig into the noise" . However, here we are limited by the small 
number of observations. Points to note are: (i) even with small n, we detect 
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Figure 2. Announced planets (open circles) and candidates (filled 
squares) (periodogram false alarm < 10~^). Sloping dashed lines show 
different values of Msini for a 1 Mq star. Vertical dashed line shows 
maximum period searched. We show only fits with Msini < WMj. 



excess variability if K > 1-2 a, (ii) it is hard to characterize an orbit when the 
number of observations is < 10, (iii) we detect periodic signals with K ~ 2-4 a 
when n = 10—20, and (iv) detecting amplitudes < 1 a requires n > 50. 

We are currently extending this approach to include eccentric orbits. How- 
ever for now, we use the best-fit circular orbit amplitude and period from the 
periodogram as initial conditions for a full Keplerian fit. The resulting Keple- 
rian amplitudes and periods for the significant detections are shown in Figure 2. 
Most, but not all, of the candidates from the periodogram analysis are announced 
planets (open circles). Those candidates that have not been announced (filled 
squares) have either low amplitude or long orbital period. There is an "effective 
detection threshold" beyond that set by the statistics of the periodogram. This 
comes from (i) we must be sure the observed variations are not due to periodic 
stellar jitter, and (ii) inability to characterize orbits longer than the duration of 
the observations. The announced planets have velocity amplitude > 10-15 m/s, 
periods less than the duration of the observations (< 3-5 years), and number of 
data points n > 10. Better characterization of jitter would aid detection of low 
amplitude planets. 
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Finally, for most stars with no detection, we can exclude companions with 
velocity amplitudes > 10-20 m/s, or Msini > (0.35-0.7)Mj {a/AUy/^, for 
orbital periods less than the duration of the observations (2 to 5 years). 

3. Constraining the Distribution of Orbital Parameters 

Wc adopt a Bayesian approach to constrain the distribution of masses, orbital 
periods, and eccentricities using both detections and upper limits. These calcu- 
lations are currently in progress, here we discuss the methodology. Consider a 
model of the distribution of masses and orbital periods n(M, P) (normalized so 
that / n{M, P)dMdP = 1) and fraction of stars with planets /. We write a 
likelihood function L = Hj Lj where j denotes each star, and 

= (1 - f)q, +fj dMdP n{M, P) pj{M, P). (1) 

This expression includes the probability of the data given a planet of mass M, 
period P, 

N 



MM.P) = I i(cos.) J- n 



d(f> YJ [ [vk - fkiMsmi,P,i 



(2) 

and the probability of the data given no planet qj. No decision must be made for 
each star about whether a planet has been detected or not: both possibilities are 
included, weighted by their relative probability. This approach can be readily 
generalized to include, for example, eccentric orbits. 

As a check, imagine we can definitely say whether a star has a planet or not. 
The likehhood is then L = /^p(i — /)^*-^p where Np planets are detected out of 
A'* stars. Maximizing L with respect to / gives the expected result / = Np/N^,, 
the best guess at the fraction of stars with planets is the observed fraction. 

Two striking features already apparent are (i) the lack of massive planets at 
short periods (easiest to detect, but not seen), and (ii) the "pile- up" of planets 
at orbital periods ~ 3 days (no planets were found between 2 and 3 days in 
our search, suggesting that this is not a selection effect). We are currently 
investigating these issues and others more carefully, including (i) the fraction of 
stars with planets, (ii) whether there is a lack of Saturn-mass companions at 
short orbital periods, and (iii) whether planet properties depend on the spectral 
type of the host star. 
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